A Model for Nonequilibrium Wetting Transitions in Two Dimensions 
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A simple two dimensional (2d) model of a phase growing on a substrate is introduced. The model 
is characterized by an adsorption rate q, and a desorption rate p. It exhibits a wetting transition 
which may be viewed as an unbinding transition of an interface from a wall. For p = 1, the model 
may be mapped onto an exactly soluble equilibrium model exhibiting complete wetting with critical 
exponents 7 = 1/3 for the diverging interface width and xq = 1 for the zero-level occupation. For 
<p / 1 a crossover to different exponents is observed which is related to a KPZ type nonlinearity. 

PACS numbers: 68.45. Gd; 05.40. +j; 05.70.Ln; 68.35.Fx 



The interaction of a bulk phase (a) of a system with a 
wall or a substrate may result in very interesting wetting 
phenomena. In particular a layer of a second phase (f3) , 
which is preferentially attracted to the wall may be 
formed in its vicinity. As some of the parameters control- 
ling the system, say temperature or chemical potential, 
are varied, the thickness of the (3 layer may diverge, lead- 
ing to a wetting transition. Such transitions have been 
theoretically studied and experimentally observed in a 
variety of systems and models in thermal equilibrium (for 
a review, see Ref. Q]). Wetting transitions may be viewed 
as the unbinding of an interface from a wall. Within this 
approach one considers the interface configuration h(r) 
which gives the height of the interface above the wall at 
point r. One then introduces an effective Hamiltonian of 
the form 0] 



n= I d^rl-iVhf + V{h{r))\ 



(1) 



where a is the surface tension of the a- (3 interface, 
V(h(r)) yields the effective interaction between the wall 
and the interface, and d — 1 is the interface dimension. 
The potential V which contains an attractive component 
may bind the interface to the wall. However, as the tem- 
perature or other parameters describing the system are 
varied, the attractive component of the potential may be- 
come weaker and it is no longer able to bind the interface, 
leading to a wetting transition. In d = 2 dimensions one 
usually distinguishes between critical wetting and com- 
plete wetting. Critical wetting is marked by the diver- 
gence of the interface width when the temperature T is 
increased towards the transition temperature TV mov- 
ing along the coexistence curve of the a and (3 phases. 
On the other hand complete wetting is characterized by 
the divergence of the interface width when the chemi- 
cal potential difference between the two phases is varied, 
moving towards the coexistence curve at T > Tw- These 
types of transitions are associated with two different sets 
of critical exponents. 



A very interesting question which has not been studied 
in detail so far is that of wetting transitions under non- 
equilibrium conditions. Here the (3 phase is adsorbed to 
the wall via a growth process whose dynamics, unlike that 
of equilibrium processes, does not obey detailed balance. 
This problem may be studied by considering the behav- 
ior of a moving interface interacting with a wall. Such 
transitions have been reported in recent studies of the 
dynamics of certain models of coupled maps || . 

In this Letter we introduce a class of nonequilibrium 
growth models of a one dimensional interface interacting 
with a substrate. The interface evolves by both adsorp- 
tion and desorption processes which in general do not 
satisfy detailed balance. By varying the relative rates 
of these processes, a transition from a binding to a non- 
binding phase is found. For a particular value of the 
desorption rate, for which the dynamics happens to have 
detailed balance, the model may be mapped onto an ex- 
actly soluble equilibrium model which exhibits a com- 
plete wetting unbinding transition. The associated crit- 
ical exponents are 7=1/3 for the interface width and 
xq = 1 for the base level occupation. For generic values of 
the desorption rate, however, detailed balance is violated 
and a crossover to different exponents is observed. 

A. Definition of the model: The model is defined 
in terms of growth of a Id interface, in which both ad- 
sorption and desorption processes take place. We con- 
sider a restricted solid-on-solid (RSOS) growth process, 
where the height differences between neighboring sites 
are restricted to take values 0, ±1. The model is defined 
on a Id lattice of N sites with associated height vari- 
ables hi = 0, 1, . . . , oo and periodic boundary conditions. 
We use random sequential dynamics which are defined 
through the following algorithm: at each update choose 
a site i at random and attempt to carry out one of the 
processes 

- adsorption of an adatom with probability q/(q+p+l): 



hi + 1 



(2) 
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FIG. 1. Phase diagram for an infinite system. The wetting 
transition takes place along the solid line. Along the dashed 
line the growth velocity in the model without a wall does not 
depend on a global tilt of the interface, indicating that the 
effective KPZ nonlinearity vanishes. Both lines intersect in 
the point q = p = 1, where the microscopic processes obey 
detailed balance. 



- desorption of an adatom from the edge of an island 
with probability l/(q + p + 1): 



h, t -> min(h i - 1 ,h i ,h i+1 ) 



(3) 



- desorption of an adatom from the interior of an island 
with probability p/(q + p + 1): 



hi - 1 if hi-i = hi = h, H 



(4) 



If the selected process would result in a violation of the 
RSOS constraint \hi — hi+i\ < 1, the attempted move is 
abandoned and a new site i is selected. In addition, a 
hard-core wall at zero height is introduced, i.e. a pro- 
cess is only carried out if the resulting interface heights 
are non-negative. One can prove that these processes in 
general do not satisfy detailed balance. 

The presence of a hard-core wall at ft, = leads to a 
phase transition that takes place even in finite systems. 
This can be seen as follows. Without the wall the in- 
terface in a finite system has a finite width. For fixed 
p > the parameter q controls the mean growth velocity 
of the interface, i.e. for large q the interface grows while 
for small q it moves downwards. These two regimes are 
separated by a critical growth rate q — q c for which the 
mean velocity is zero. Therefore, on large time scales, a 
lower wall will only affect the interface dynamics if the 
interface does not move away from the wall, i.e. q < q c , 
resulting in a smooth interface. In the growing phase 
q > q c , however, the interface does not feel the wall. It 
is rough and propagates with a constant mean velocity. 
The phase transition line for an infinite system is shown 
in Fig. 1. 



Throughout this paper we are particularly interested 
in the the mean growth velocity v in the growing phase, 
the occupation po of the zero-height layer in the smooth 
phase, and the interface width in the smooth phase which 
is defined by 



w 



N N 

±y( hi -±y h 



(5) 



Consider now the thermodynamic limit N — ► oo. Near 
criticality in the growing phase q > q c , we expect the 
interface velocity to scale like 



« ~ (Q ~ <lc) v , 



(6) 



whereas in the smooth phase, q < q c , the expected scaling 
for bottom layer occupation and width is 



Po ~ (q c - q) 



■> - n 



(?c - q)~ 



(7) 



It would be interesting to find out how the critical ex- 
ponents y, X(s, and 7 depend on p. We note that the 
case p = is special: In this case atoms cannot be des- 
orbed from a completed layer, and the interface cannot 
move below its actual minimum height. This means that 
the hard-core wall becomes irrelevant. Therefore, the 
phase transition (which still exists for p — 0) relies on a 
completely different mechanism. The p = transition is 
expected to belong to the universality class of a closely re- 
lated model previously considered in Ref. H . It has been 
shown that in this case some of the critical exponents can 
be related to the universality class of directed percola- 
tion (DP). In particular one expects xq = and y = v±, 
where f3 = 0.276 and v± = 1.73 are the density and 
correlation length exponents, respectively, of DP. Here 
completed layers play the role of absorbing states from 
where the system cannot escape. For p > 0, however, the 
system is ergodic and one expects different critical expo- 
nents. Another special case is p = 1: here the system 
does satisfy detailed balance and can be solved exactly. 
In the following, we present our analysis of the wetting 
transition for p = 1. We then consider the general case 
for which the model does not have detailed balance. 

B. Exactly soluble case, p — 1: We first show that in 
this case the steady state satisfies detailed balance and 
that for q < 1 the probability of finding the interface in 
a particular configuration an = {hi, . . . , ftjv} is given by 
the distribution 

P(h u ...,h N ) = P(a H ) = Z N X q H{^,..,h N ) i (8) 
where 



N 



H = H(h 1 ,...,h N ) = Y,h i 



(9) 



i=l 



is the sum of all heights. Here, the partition sum Zn — 
Y] h Hn q H runs over all interface configurations which 
respect the RSOS and the hard-core wall constraints. 
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In order to prove Eq. (||), notice that the processes 
(||)-(|J) subjected to the RSOS constraint correspond to 
a change of H by one unit. Therefore, if the distribution 
of Eq. (||) is to hold in the steady state, the probabilities 
of finding the interface in states with total height H and 
H + 1 have to satisfy 



P(a H+1 )/P(a H ) = q. 



(10) 



As one can read off from the processes (||)-(||), any al- 
lowed transition <jh - ► o~h+i occurs with rate w(<7h —> 
<Jh+i) = q- Moreover, for each such process an —* o-r+i, 
there is a reverse process <jh+i ~~ * °~h which takes place 
with a nonzero rate. For p — 1, this rate is given by 
w(uh+i —* o"ff) = 1 so that 

w{cr H —> a H+1 )/w{a H +i -> <Jh) = q- (p = 1) (11) 

Together with Eq. (|Io| ) this implies that the processes 
(||)-((|) satisfy detailed balance. Notice that the above 
consideration is consistent with the hard wall constraint 
hi > 0. Equation (||) yields the steady state distribution 
only for q < 1. The unbinding transition takes place at 
q = 1 and the height distribution becomes time depen- 
dent for q > 1. For p ^ 1, detailed balance is violated. 
This can be proven by construction of explicit cycles of 
configurations in small systems for which the rates of 
moving clockwise and counter-clockwise are unequal. 

We now apply the transfer matrix formalism m to 
study the distribution (||). Let us define a transfer ma- 
trix T acting in spatial direction by 



Th,h' — 



q" 




i£\h-h'\ < 1 
otherwise 



(12) 



where h, h' > 0. Steady state properties can be derived 
from the eigenvector 4> that corresponds to the largest 
eigenvalue p of the transfer matrix Y^h'=o Th,h'4>h' 
ji (j) h . From the squares of the eigenvector components 
one can derive various steady state quantities. For exam- 
ple, the probability ph of finding the interface at height 
h is given by pn — 4>h/ J2h' ^h' ■ Here we are particu- 
larly interested in the scaling behavior of bottom layer 
occupation po and the width w 2 — ^2 h {h — h) 2 where 
h = hph denotes the mean height. 

Close to criticality, where e = 1 — q is small, one can 
carry out the continuum limit 4>h — * 4>{h), replacing the 
discrete heights h by real-valued heights h. Then, the 
above eigenvalue problem turns into a differential equa- 
tion H which, to leading order in e, is given by 



(J^ + (3- M )-3efc) cf>(h) = 0. 



(13) 



This equation, together with the boundary conditions 
4>(0) = <H°°) = 0; has a unique physical solution. Sim- 
ple dimensional analysis indicates that the height vari- 
ables scale as h ~ e -1 / 3 and thus the width diverges 
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FIG. 2. Results obtained from Monte Carlo simulation of 
the growth model with 1500 sites. The width w and the bot- 
tom layer occupation po are measured for (*) p=0, (o) p=0.05, 
(•) p=l, and (o) p=2.0. For p=l we fitted straight lines. The 
curvature for p / 1 indicates a crossover to KPZ exponents. 



as w 2 ~ e -2 / 3 . The occupation of the bottom layer 
in the continuum limit is given by po = A/" _1 (0'(O)) 2 , 
where J\f = f dh(f> 2 (h) is a normalization factor. Since 
</>'(0) ~ e 1 / 3 and N ~ e" 1/3 one obtains a linear scaling 
law po ~ e. The critical exponents for p — 1 are thus 
given by 

x = l, 7-1/3- (14) 

C. Numerical results: In order to determine the crit- 
ical exponents for other values of the growth rate p, we 
perform Monte-Carlo simulations. The width w and the 
occupation of the bottom layer po are measured in the 
smooth phase q < q c . Depending on e = q c — q, we first 
equilibrate a system of size 1500 over a time interval up 
to 4 x 10 6 time steps. Then the thermal averages of w 
and po are measured over a time interval of the same 
size. Similarly, the interface velocity v is measured in 
the growing phase q > q c . 

The numerical data measured in the smooth phase are 
shown in Fig. 2. From the slopes in the double loga- 
rithmic plots we estimate the critical exponents (see Ta- 
ble I). For p = 1, the numerical results we obtain are 
consistent with the exact values derived above. In ad- 
dition the velocity exponent is found to be y = 1.01(3). 
For < p < 1 we observe a crossover to different critical 
exponents. Since the values differ from those in Eq. (|l4| ) 
by less than 30% and the crossover is extremely slow, it 
is difficult to determine these exponents precisely. Our 
estimates are y — 1.00(3), xo — 1-5(1), and 7 = 0.41(3). 
Similar results are obtained for p > 1 except for xo which 
is close to one in this case. Finally, for p — ► 0, we observe 
another crossover. This is consistent with the results of 
Ref. which indicate that the interface width diverges 
logarithmically at the p = transition. 

D. General Considerations: Equilibrium wetting of a 
2d Ising system has been studied in Refs. It was 
shown that a column of weak bonds (which acts as an 
attractive potential for the domain wall separating the 
up and down states) located at the boundary induces 
a wetting transition at some finite temperature Tw- At 
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v = 


p = 0.05 


p = 1.0 


p = 2.0 


Qc 


0.3991(5) 


0.5564(2) 


0.9999(2) 


1.2329(5) 


x 


0.27(1) 


1.51(6) 


0.96(5) 


1.02(5) 


7 


(log.) 


0.41(3) 


0.32(3) 


0.37(3) 


y 


1.69(5) 


0.98(3) 


1.01(3) 


1.00(3) 



TABLE I. Estimates for the critical exponents. 



the critical wetting transition the interface width diverges 
with 7 = 1. 

The transition found for p = 1 is of different nature. 
In this case, (1 — q) acts as a chemical potential differ- 
ence between the two coexisting phases. For q < 1 the 
chemical potential difference drives the interface towards 
the wall, resulting in a smooth interface. On the other 
hand, for q > 1 the interface is driven away from the wall, 
resulting in a KPZ-like rough moving interface J§|. The 
critical behavior associated with the transition is thus 
that of complete wetting. 

For p ^ 1, the mapping to equilibrium is impossible 
since detailed balance is violated. Here a KPZ-type non- 
linearity is expected to be responsible for the different 
exponents we observe. Within this approach, one de- 
scribes the system by the Langevin equation 

dh(r,t) „ 2l . , dV(h(r,t)) 

+A(VMr,i)) 2 +CM), 

where £( r j0 is a zero-average Gaussian noise field with 
variance (C(r, t)C(r J , t')) = 2DS d - 1 (r - r')5(t - t') and V 
is the effective interaction between the wall and the inter- 
face. This equation has been studied recently ||[l0) in the 
context of nonlinear diffusion with multiplicative noise. 
A simple scaling argument suggests that the width expo- 
nent corresponding to the wetting transition described 
by this equation is given by 

7 = (2-z)/(2z-2), (16) 

where z is the dynamic exponent. For a Id interface 
z = 3/2, yielding 7 = 1/2. Our numerical results in- 
dicate that the width exponent 7 is larger for p =/= 1 as 
compared with its 1/3 value at p = 1, although it seems 
to be smaller than 1/2. However, in view of the very slow 
crossover expected in this problem (see below) it is possi- 
ble that 7 is indeed 1/2, but more extensive simulations 
close to q c would be needed to observe it. 

The bottom layer occupation po may be related to £ , 
where £ is the correlation length. Thus, we expect xq to 
be equal to the correlation length exponent u, which for 
the KPZ equation is given by v = l/(2z — 2), yielding 
xq = 1 in Id. However, the numerical results suggest that 
this scaling argument is valid only for p > 1 , whereas for 
< p < 1 much larger values for xq are obtained. This 
may be related to the existence of different universality 



classes in both cases, corresponding to the distinction 
between an 'upper' and a 'lower' wall in Ref. ppf . 

In order to verify this picture we made a numerical es- 
timate of the effective nonlinear KPZ term corresponding 
to the RSOS model considered in this work. This is done 
by comparing the growth velocities of a flat and a tilted 
interface in absence of a wall. We find that the nonlincar- 
ity is indeed non vanishing in the (p, q) plane, except on 
a particular line (the dashed line in Fig. 1). As expected, 
this line and the phase transition line are different and in- 
tersect in the point q = p = 1. At all other points on the 
transition line the nonlinear term is not vanishing, and 
the KPZ-like exponents are expected to be valid. Since 
A is small in the vicinity of the q = p = 1 point, very 
slow crossover phenomena occur, making it difficult to 
observe the true exponents in this region. Moreover, the 
interaction between the wall and the interface strongly 
depends on the sign of A. This may lead to to different 
exponents xq for the bottom layer occupation on the two 
sides of the special point. 
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